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The density-matrix renormalization group is used to study the phase diagram of the one- 
dimensional half-filled Hubbard model with on-site (U) and nearest-neighbor (V) repulsion, and 
hopping t. A critical line V C (U) Ri U/2 separates a Mott insulating phase from a charge- density- 
wave phase. The formation of bound charge excitations for V > 2t changes the phase transition from 
continuous to first order at a tricritical point Ut ~ 3.7t, Vt = 2t. A frustrating effective antiferro- 
magnetic spin coupling induces a bond-order-wave phase on the critical line V C (U) for Ut < U < It. 
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The properties of quasi-one-dimensional materials have 
been extensively studied in recent years El @, [|. These 
materials exhibit rich phase diagrams and display un- 
usual optical properties due to the combination of re- 
duced dimensionality and strong electronic correlations. 
Consequently, much effort has been devoted to under- 
standing the ground-state and optical properties of theo- 
retical one-dimensional correlated electron systems such 
as the half-filled Hubbard model with on-site (U) and 
nearest-neighbor (V) repulsion and hopping term t. Nev- 
ertheless, the ground state phase diagram of this model 
is still controversial §, |, |, @, §, §, |§. It is known 
H that the system is a Mott insulator for U > 2V and 
a charge density wave (CDW) insulator for U < 2V. 
The quantum phase transition is continuous at weak 
coupling (U, V <C t) and first order at strong coupling 
(U, V » t). Numerically [§, |, 0, one finds that the or- 
der of the transition changes at a tricritical point (Ut, Vt) 
with Vt/t « 1.5 — 2.5, but this feature is not well under- 
stood. Recently, it has been proposed || ^, |l(| that a 
bond-order-wave (BOW) phase exists between the Mott 
and CDW phases up to the tricritical point. Concur- 
rently, the optical properties of this model have been de- 
termined in the Mott insulating phase 0, [l^, 0, f[5| . 
In particular, it has been found that the lowest optical 
excitations consist of a pair of independent charge excita- 
tions for V < 2t, while they are bound states for V > 2t. 
Surprisingly, the remarkable proximity of the tricritical 
point to the boundary between bound and free charge 
excitations has not been noticed until now. 

Here I investigate the ground-state phase diagram of 
the half-filled one-dimensional extended Hubbard model 
in the repulsive regime using the density-matrix renor- 
malization group (DMRG) [^6|. I show that the nature 
of the low-lying charge excitations determines the order 
of the transition and the position of the tricritical point. 
A BOW phase is found only at intermediate coupling on 
the critical line V C (U) between Mott and CDW phases. 

The model is defined by the Hamiltonian 



+V£(ft 1 -l)(ft, +1 -l). (1) 
I 

Here ct Q.o- are creation and annihilation operators for 
electrons with spin a =],[ at site I — 1, . . . , N, hi >a = 
ct^c^cr, and hi = nj f + n;^. I exclusively consider sys- 
tems with an even number N of sites. At half filling the 
number of electrons equals N. The interaction is repul- 
sive U > V > 0. The Hamiltonian ([j]) has a particle-hole 
symmetry and a spatial reflection symmetry. Thus, each 
eigenstate has a well-defined parity under charge conju- 
gation (P c = ±1) and belongs to one of the two irre- 
ducible representations, A g or B u , of the reflection sym- 
metry group. The ground state belongs to the symmetry 
subspace = (A g ,P c ) and optically excited states be- 
long to the symmetry subspace B~ = (B u , —P c ) because 
the current operator is antisymmetric under charge con- 
jugation and spatial reflection |l4[ . 

DMRG Jl6| is known to be a very accurate numeri- 
cal method for one-dimensional quantum systems with 
short-range interactions such as the Hamiltonian (jl]). 
Here I use the finite-system DMRG algorithm to cal- 
culate eigenenergies and static properties at low-energy. 
Spin gaps E s = E (N, 1) — E o (N,0) and single-particle 
charge gaps E c = 2[E (N + 1, 1/2) - E (N,0)} are de- 
rived from the ground state energies E (N e ,S z ) for a 
given number of electron N e and a given spin S z . I also 
use a symmetrized DMRG Jl5|, method to calculate 
the energy and static properties of the lowest eigenstates 
in the B~ symmetry subspace. 

All DMRG methods have a truncation error which is 
reduced by increasing the number m of density-matrix 
eigenstates kept in the renormalization procedure [ fl6|| . 
To achieve a greater accuracy and to obtain error esti- 
mates I extrapolate DMRG results to the limit of van- 
ishing discarded weight P m . For eigenenergies DMRG 
errors vanish linearly with P m , while for other quantities 
DMRG errors usually scale as (P m ) a with < a < 1. 
The largest number of density-matrix eigenstates used in 
this work is m = 1200 and truncation errors are negli- 
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TABLE I: Critical nearest-neighbor repulsion V C (U), single- 
particle charge gap E c at (U, V c ), and transition type for sev- 



eral values of U. 


Numbers in parenthesis are 


estimated errors. 


U/t 


Vc(U)/t 


E c /t 


transition order 


2 


1.125 (25) 





continuous 


3 


1.640 (10) 





continuous 


4 


2.150 (10) 


0.05 


first order 


5 


2.665 (5) 


0.11 


first order 


6 


3.155 (5) 


0.15 


first order 


8 


4.141 (4) 


0.45 


first order 


12 


6.115 (5) 


2.0 


first order 


40 


20.041 (4) 


29.3 


first order 



gible for all results presented here. DMRG calculations 
have always been carried out for several sizes N in order 
to check finite-size effects and, if necessary, results have 
been extrapolated to the thermodynamic limit N — ► oo. 
The largest system size used in this work is N = 1024. 

For V < U/2 the model (ffl) describes a Mott insulator. 
[This phase is often called a spin-density-wave (SDW) 
phase because antiferromagnetic spin correlations decay 
algebraically in the ground state.] There is no broken 
symmetry and the ground state is nondegenerate. The 
charge gap E c is finite but the spin gap E s vanishes in 
the thermodynamic limit. In this Mott insulating phase 
the lowest eigenstate in the B~ symmetry sector always 
contributes to the optical spectrum. Therefore, the dif- 
ference between the B~ and A+ ground-state energies 
corresponds to the optical gap E opt . For V > U/2 the 
system is in a long-range ordered CDW phase with a 
doubly degenerate ground state in the thermodynamic 
limit. In this phase both charge and spin gaps are fi- 
nite. The CDW order parameter < \m e \ < 1 gives 
the amplitude of the ground-state density modulation 
(hi) = 1 + (— l)'m e . Such a CDW breaks the charge- 
conjugation symmetry and the spatial reflection symme- 
try. Thus one of the two degenerate ground states be- 
longs to the A+ symmetry subspace and the other one 
to the B~ subspace. Obviously, this consideration shows 
that the low-energy B~ excitations in the Mott phase 
play an important role in the Mott-CDW transition. 

The transition from the CDW phase to the Mott phase 
can be investigated with DMRG as done previously for 
the CDW-metal transition in the Holstein model Jl8| . 
More precisely, I have investigated the lowest eigenstates 
to determine the ground-state degeneracy and symme- 
try and I have calculated the electronic staggered sus- 
ceptibility to check if the ground state has long-range 
CDW order. Some results for the critical line V C (U) 
determined with this approach are given in Table §. 
These results agree quantitatively with recent quantum 
Monte Carlo (QMC) simulations and the Gaussian 
transition line in Ref. ||, although the interpretation of 
this phase boundary is completely different (see below). 
These values are also in good agreement with strong- 
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FIG. 1: Ground-state double occupancy d versus V close to 
the critical coupling V C (U) for several U (in units of t). 

coupling perturbation theory || down to U = 6t and the 
results of early numerical simulations |Q, This con- 
firms the accuracy of the present approach and I will not 
elaborate further on the location of this boundary. 

To determine the order of the transition on the crit- 
ical line V C (U), one can examine the derivatives of the 
ground-state energy per site with respect to the interac- 
tion parameters U and V. Using the Hellmann-Feynman 
theorem one shows that the derivatives are given (up 
to known constants) by ground-state expectation values 
of the double occupancy and nearest-neighbor density- 
density operators, respectively. These expectation val- 
ues can be calculated accurately with DMRG. Figure [j] 
shows the average ground-state double occupancy 

d = jr^2{nnhn) (2) 
I 

versus V close to the critical coupling V C (U). For U — 
404, d jumps from almost zero to slightly less than 0.5 at 
V c and thus the transition is clearly first order, while it 
is continuous for U — 2t. For U — 6t the derivative of 
d with respect to V evolves into a (5-function at V = V c 
as the system size N increases. Therefore, d is discon- 
tinuous at V C (U) and the transition is also first order for 
U = 6t. With this method I have found that the transi- 
tion at V C (U) is first order from the strong-coupling limit 
(U, V > t) down to at least U = At (V c ~ 2.154) and is 
continuous from the weak-coupling limit (U, V -C t) up 
to at least U = 3.5t (V c w 1.9*). 

An investigation of the low-lying charge excitations in 
the Mott insulating phase allows us to understand the 
existence of a tricritical point and to determine its posi- 
tion. Figure || shows the evolution of the charge gap E c 
and the optical gap -Eopt for three different couplings U. 
I have always found that the charge gap E c varies con- 
tinuously and goes through a minimum at V C (U) as V is 
increased for fixed U. For weak coupling (U < 3t) the gap 
vanishes smoothly on the critical line. For stronger cou- 
pling (U > 4) the slope of E c with respect to V becomes 
very large and is discontinuous at V c . In this regime the 
charge gap remains finite at the transition (see Table |) . 
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FIG. 2: Charge gap E c (upper line) and optical gap E opt 
(lower line) versus V close to the critical coupling V C (U) for 
U/t = 40 (solid), 6 (dashed), and 2 (dot-dashed). For U = 2t, 
E op t = E c . E op t is not shown in the CDW phase (V > V c ). 



In a Mott insulator, elementary charge excitations can 
be understood as spinless bosons in Hubbard bands. Ex- 
cited states in the B~ subspace always consist of an even 
number of such elementary excitations with zero total 
charge [[ll| |l3[ DMRG calculations combined with 
analytical methods |l~4|, |l5| reveal their properties. For 
V < 2t the low-energy B~ excited states consist of two 
independent elementary charge excitations and the op- 
tical gap equals the charge gap in the thermodynamic 
limit. For V > 2t, however, the lowest B~ excitations are 
bound states starting at an energy E opt < E c (see Fig. ||). 
Just above 2t they are excitons made of two elementary 
charge excitations. Close to the critical line V C (U) ~ U/2 
they become finite-size "droplets" of the CDW phase, 
(i.e., a bound state of several elementary charge excita- 
tions). The excitation energy E opt of the lowest CDW 
droplet remains finite at the transition V C (U) at least for 
U > 5 (see also Fig. §. CDW droplet sizes increase 
sharply as V approaches the critical coupling V C (U) for 
fixed U and reach a finite value £ c at the critical line. This 
critical droplet size £ c diverges for U — > oo but tends to 
2 (corresponding to an exciton) if V c approaches 2t. 

Combining all results, one discerns two distinct regimes 
in the Mott phase as the critical line V C (U) is approached. 
If V c < 2t, the transition is continuous, the B~ exci- 
tations are made of free charge excitations and become 
gapless at V c . If V c > 2t, the transition is first order, the 
low-lying B~ excitations are bound states and remain 
gapped at V c . Therefore, I believe that the tricritical 
point is precisely located at the intersection of the V = 2t 
boundary with the quantum critical line V C (U) ~ U/2, 
i.e., at U t ~ 3.7i,Vt = 2t. 

Adapting Hirsch's analysis 0, one can now qualita- 
tively explain why the transition from the Mott to the 
CDW phase changes from first order to continuous. If 
V C (U) > 2t, CDW droplets with sizes £ > £ c are energet- 
ically favored when V becomes larger than V C (U), and 
the system tunnels to the CDW phase by nucleation. On 
the critical line V C (U) the creation energy E opt and crit- 



ical size £ c of the lowest CDW droplet diminishes with 
U. At V c = 2t, CDW droplets (with sizes £ -> 2) become 
instable and split into independent gapless elementary 
charge excitations. For V C (U) < 2t, it is then advanta- 
geous to create many such excitations as V is increased 
above V C (U) and the transition becomes continuous. 

In the field theory approach to the lattice model (|l|) 
the limit between free and bound charge excitations cor- 
responds to a parameter K p = 1 /2 for a half-filled system 
p3| , p"4| ]. Field-theoretical predictions for the low-energy 
optical spectrum agree perfectly with dynamical DMRG 
calculations |l4|, 15 . Thus the tricritical point probably 
corresponds to K p — 1/2. In Ref. 10, it has also been 



found that K p is close to 1/2 in this region. 

Recently, Nakamura |§ has proposed that a narrow 
region with long-range BOW order exists between the 
CDW phase and the Mott phase at weak to intermediate 
coupling. The BOW region vanishes at the tricritical 
point, beyond which there is a direct first-order transition 
from the Mott phase to the CDW phase. In this scenario 
the critical line V C (U) calculated above would correspond 
to the BOW-CDW phase boundary for U < U t . The 
existence of long-range BOW order in the critical region 
has been confirmed by QMC simulations |[oj . 

In a BOW phase the ground state is two-fold degener- 
ate and both charge and spin gaps are finite. A long- 
range ordered BOW corresponds to a finite staggered 
bond order 

\ E ( d t d '+^ + e m,<A-) = po + (-!)'<* (3) 

with 5 =/= in the thermodynamic limit. To distinguish 
the BOW phase from the Mott phase I have calculated 
the spin gap E s and the bond order parameter 5 using 
DMRG. After extrapolation to vanishing truncation er- 
rors P m and to the thermodynamic limit N — > oo, this 
approach gives spin gaps with an accuracy of 10~ 3 t or 
better and allows one to detect a bond order alternation 
as small as \5\ = 0.01. (It appears numerically that 5 
either vanishes as 1/yN or tends to a finite value with 
finite-size corrections scaling as 1/N.) I have found a 
BOW ground state only in a very narrow region adja- 
cent to the critical line V C (U) for intermediate coupling 
At < U < 6t. For U = At and V = 2.1At, I have ob- 
tained S = 0.08 in quantitative agreement with QMC 
simulations ]l0|j . However, the BOW phase extends to 
significantly stronger coupling than reported in Ref [lC| . 
For instance, there is a BOW ground state with 5 « 0.12 
for U = 6t and V — 3.145i. I have not observed any 
BOW ground state at weak coupling (U < 3i), at strong 
coupling (U > 8t), and for V < U/2. Actually, as BOW 
ground states are found only within 0.02t of the criti- 
cal line V C (U), I believe that the BOW phase exists only 
on this critical line for intermediate coupling U starting 
from the tricritical point Ut ~ 3.7£ up to a upper limit 
which is smaller than 8t. Field theory H] also suggests 
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FIG. 3: Schematic ground-state phase diagram: the tran- 
sition is continuous at weak coupling (dashed line) and first 
order at strong coupling (both solid lines). A circle marks 
the tricritical point. The dot-dashed line is the boundary be- 
tween free and bound charge excitations in the Mott phase. 
The thick solid line indicates the BOW phase and a square 
marks the transition from a "spin-liquid" to a dimerized spin 
ground state on the critical line. 



suits of most previous investigations of this problem (see 
Refs. |] to [?] and references therein). 

In summary, I have studied the ground-state phase di- 
agram of the one-dimensional half-filled extended Hub- 
bard model for U, V > using DMRG. I have shown that 
the phase diagram is determined by 3 mechanisms: First, 
the competition between the on-site (U) and nearest- 
neighbor (V) repulsion is responsible for a quantum crit- 
ical transition from a Mott to a CDW ground state at 
V ~ U/2. Second, the competition between the hopping 
term t and the nearest-neighbor repulsion V is responsi- 
ble for the transformation of free charge excitations into 
bound states (CDW droplets) at V = 2t in the Mott 
phase and thus changes the Mott-CDW transition from 
continuous to first order. Third, a frustrating effective 
antiferromagnetic spin coupling leads to a dimerized spin 
ground state (and thus to a BOW phase) on the Mott- 
CDW critical line at intermediate coupling. 



that a BOW phase can occur only on the boundary be- 
tween CDW and Mott phases in the model (H). Even in 
the Hartree-Fock approximation, one finds a BOW phase 
only on the critical line U — 2V between a SDW phase 
and a CDW phase Jl|. 

The occurrence of the BOW phase on the critical line 
for U < 8t can be understood as the result of increasing 
frustration in the spin degrees of freedom. On this criti- 
cal line, low-energy charge excitations are dispersionless 
and thus ineffective in a degenerate CDW-Mott ground 
state. A strong-coupling expansion up to 4th order in 
t/U shows that the spin properties are determined by an 
effective Heisenberg Hamiltonian with nearest- (J\) and 
next-nearest-neighbor (J2) antiferromagnetic exchange 
couplings fl. This strong-coupling perturbation anal- 
ysis gives accurate results for the critical line V C (U) (and 
thus is expected to be valid) down to U = 6t. The ra- 
tio J2/J1 is strongly enhanced by the nearest-neighbor 
repulsion V. It is known (2(], ||l[ that the frustration 
due to the J2 coupling causes a quantum phase transi- 
tion from a "spin-liquid" ground state for 4 J2 < J\ to a 
dimerized spin ground state for 4 Ji > J\ in the Heisen- 
berg model. On the critical line V C (U) of the extended 
Hubbard model, this corresponds to the appearance of a 
BOW phase (driven by the spin dimerization) for U < It, 
in agreement with my DMRG results, if one uses the val- 
ues of J\ and J2 given by the 4th-order perturbation 
expansion 0. Away from the critical line [V < V c (£/)], 
one observes only algebraically decreasing BOW corre- 
lations because charge fluctuations, which are no longer 
dispersionless, destroy the local spin moments and thus 
prevent any long-range spin order. 

Figure || shows the phase diagram found with DMRG. 
This phase diagram disproves recent speculations about 
an extended BOW phase from weak coupling up to the 
tricritical point || || [l(| but is compatible with the re- 
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